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The electronic structure in the vicinity ol the 1-heptagonal and 1-pentagonal defects in the carbon 
| graphene plane is investigated. Using a continuum gauge field-theory model the local density of 

^s^j . states around the Fermi energy is calculated for both cases. In this model, the disclination is 

represented by an SO (2) gauge vortex and corresponding metric follows from the elasticity properties 
of the graphene membrane. To enhance the interval of energies, a self-consistent perturbation scheme 
is used. The Landau states are investigated and compared with the predicted values. 
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I. INTRODUCTION 



Nanostructured carbon materials are the materials with a special geometrical structure of their molecules which 
we call carbon nanoparticles. This geometrical structure is accompanied by topological defects in a hexagonal plane 
lattice called graphene. In most cases, these defects originate from the presence of the pentagons for the positive 
, curvature and the heptagons for the negative curvature [l|. 

There are known variously shaped carbon nanostructures. The most famous is the fullerene which has the 
structure of the soccer ball and can be approximated as a sphere. It is composed of 60 carbon atoms which 
create 20 hexagons and 12 pentagons @, However, other structures exist, for example, nanocones, nanotoroids, 
nanotubes, nanohorns etc. Wide varieties of electronic properties of these structures was studied. Some examples are 
given by magnetic properties 0, HI , optical absorption spectra or electronic properties of nanotube caps @. These 
properties give a potential use in nanoscale devices like quantum wires, nonlinear electronic elements, transistors, 
molecular memory devices or electron field emitters. From the theoretical point of view, it was predicted and 
experimentally verified by scanning tunelling microscopy that metallic or semiconducting properties of carbon nan- 
otubes depend on whether or not the difference n—m of the components of the chiral vector (n, m) is multiple of 3 ■ 



More complicated structures can arise when two variously shaped parts of nanoparticles are connected by a 
region with pentagon-heptagon pairs. It seems that the best approximation for pentagonal and heptagonal areas is 
hyperboloid - positively curved for pentagons and negatively curved for heptagons [Toj . 

The electronic properties of these structures can be explored by solving Dirac equation at a curved surface 



In this paper, after introducing the computational formalism, some geometrical properties of the defects 
are researched. After doing this, we research the Gaussian curvature of the surface. Then with the help of the 
Dirac equation the local density of states (LDoS) for disclinated areas near the Fermi level (close to the zero 
energy) is calculated for both pentagonal and heptagonal defects. It will be influenced by a uniform magnetic 
field. Then we compare the electronic properties of both models and, finally, we research the corresponding Landau 
states and compare them with the approximation formulas from the earlier works. The model described in [10( is used. 
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II. BASIC FORMALISM 



First we introduce the Dirac equation in (2+1) dimensions. It has the form 

il a e» a [V^-i atl -iA^ = E^. (1) 
The wave function tp, so-called bispinor, is composed of two parts, 

ipA 



each corresponding to different sublattices of the curved graphene sheet. The gauge field a^, fj, = £,ip arises from spin 



rotation invariance for atoms of different sublattices A and B in the Brillouin zone 12 1. 

The zweibein e a stands for incorporating fermions on the curved 2D surface and it has to yield the same values of 
observed quantities for different choices related by the local SO (2) rotations, 

e a ^e' a =A^, AfeSO(2). (3) 

For this purpose, a covariantly constant local gauge field w M is incorporated jl3l |: 

M " IV« + {u>J%4 = 0, (4) 

where V ^ is the Levi-Civita connection coming from the metric g^ v (see below). Then w M is called the spin connection. 
Next, the covariant derivative is defined as 

V„ = 0„ + %, (5) 

where 

^ = \uf[ la , le ] (6) 

denotes the spin connection in the spinor representation. The Dirac matrices "f a can be replaced in two dimensions 
by the Pauli matrices a a - 

7i = -"2, 72=0-1- (7) 

is the vector potential arising from the external magnetic field. 

The metric g M „ of the 2D surface comes from the parametrization with the help of two parameters £, (p: 

(£, <p) -> t = (z(£, tp), y(£, if), z(£, <p)), (8) 

where 

< £ < oo, < tp < 2ir. (9) 

The 4 components of the metric are defined as 

9„u = d^djt. (10) 

The above-mentioned hyperboloid geometry has for both cases, heptagons and pentagons, very similar but not the 
same parametrization. We research it in separate chapters. The non-diagonal components of the metric are 

9fy = 9<pt = 0. (11) 
For the zweibeins and the diagonal components of metric holds 

e \ = tp, el = -y/g^smip, (12) 

e f = VdliSinip, el = ^/g^costp, (13) 



3 



and for the spin connection coefficients w p 



de 1 = —uj 12 A e 2 , de 2 = —uj 21 A e 1 , 



uj 12 = -uj 2 \ 



so 



Then the coefficients VL^ are 



Substituting 



and making up ([T]) we get 



12 21 i °$y/9<p<p 



/5« 



wf = uf = 0. 



f2^ = 0, fl v = iuxjz. 



1 ( S(0«f j , ; _ () , , 



d^u - (j + 1/2 - a v + A v ) J— u = 

6W 



- (j + 1/2 - + A v ) x j^-v = E^/g^u. 

9w 



(14) 

(15) 
(16) 

(17) 
(18) 

(19) 
(20) 



III. GEOMETRICAL PROPERTIES 



To find the solution of (|19|). (|20|) . the knowledge of the influence of the defects on the components of the metric 
g^ is needed. It is characterized by the Frank index v which depends on the number of the defects. 

It is possible to try to approximate the geometry by the metric of the cone But this procedure causes dis- 
crepancy in the calculation of the corresponding gauge flux [l(| ■ Here is suggested a method, how to avoid this problem. 



A. Heptagonal defects 



In case of the negative curvature and affiliated heptagonal defects, in ([5]) the parametrization seems to be 

{£,(p) — > (a cosh £ cos i^, acosh^sint^, csinh£), (21) 
where a and c are some dimensionless parameters. The corresponding diagonal components of the metric are 

o 2 2 2 2 2 / \ 

g^ = a sinh £ + c cosh £, g vv = a cosh £ (22) 

and the nonzero spin connection term 

a sinh £ 



-i 2 = l 



(23) 



The defect arises by the so-called cut and glue procedure - we cut a line in the graphene plane, add a 60° area and 
glue the arising borders 9j . The geometrical properties of the new surface can be described with the help of the gauge 

potentials which change the initial components of the metric (now denoted g$) [Uj]: 



9$ -> 9/u> =V fl K(o) 



(24) 
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where 

V M it (0) = d^ {0) + , t (0) ] . (25) 

Then 

= 9^(0) ■ du~3(o) + dp^(o)\V^v !^(o)]+ 

+ ^ ( o,[<U<o>] + (^V^o) - (<%))(^(o)) (26) 
and the components of the metric and the spin connection term will be changed in the way that 

= a 2 sinh 2 £ + c 2 cosh 2 £, g vv = a 2 a 2 cosh 2 f , (27) 

n aa sinh £ 

w« = l ^=A, a=l + i/, (28) 



where v — N/6 is called the Frank index and N is the number of heptagons in the defect. In this paper, we take 
N = 1. Let us stress that for the higher number of defects, the geometrical structure is more complicated and we 
have to take into account next assumptions 0, [l5| . 

We can encircle the origin of the defect (£ = 0) by a closed loop C c and integrate over it. The result is 

d~£ = 2-kv. (29) 



No transformation of variables can change this value. If the values of the gauge field ^ are 

W Wi=i,2 =0) H ^ )i=3 = w K0) ) (30) 

where 



= -vy/r 2 , W^ 0) = vx/r 2 , r = y / x 2 +y 2 , (31) 

then 

I d~f = 2-kv= I W^dx", (32) 

so ^ serves as a vortex-like potential with a nonzero flux. This flux should be eliminated by the corresponding 
integral over the spin connection, so we must get 

lim j> u 12 d<p = -2-kv. (33) 

Substituting (1281) into the appropriate integral, the required result is really obtained. 

For our purpose, the gauge field a v — N/4. In the general case, a v depends on two constants, N and M as 
a v = N/4+M /3, where M — — 1, 0, 1 for an even number of defects and M = for an odd number of defects [9l. Hoi. fl2j . 

If the magnetic field is chosen in such a way that A = B(y, —x, 0)/2, then 

A v = -$cosh 2 £, A c = 0, (34) 

where 

$ = \a 2 <s> B, $ = (35) 

I nc 



Geometric units are used, i.e. e — h — c = 1. 
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B. Pentagonal defects 

The case of the positive curvature is described in more detail in [Toj . The parametrization is changed into 

(£, <f) — > (a sinh£ cos ip, a sinh£ shi(y9, ccosh£), (36) 
and the diagonal components of the metric 

2 2 2 2 22 / \ 

g^ = a cosh £ + c sinh £, g vv = a sinh £. (37) 

— k(0) 

Introducing the gauge potentials as for the heptagonal defects, the component g vv of the metric changes in 

the way that 

g v <p = a 2 a 2 sinh 2 £, (38) 

where a = 1 — v. This means that in the cut and glue procedure we cut a 60° area instead of insertion. Then the 
nonzero spin connection term is 

u p = 1 _ oocoshC (;j9) 



The values of the gauge field and the magnetic field are the same as in the previous case, 

a v =N/A, ~A = B(y, -x,0)/2, (40) 

so for the parametrization chosen 

A v = -$sinh 2 £, A e = 0, (41) 

where $ is defined as in (|35|) . 

IV. CURVATURE 

The Gaussian curvature is defined as 

K = (d xx f)(dyyf) ~ (0 X yf) 2 

(1 + (d x f) 2 + {dyf) 2 ) 2 ' 1 ) 

where / means the z coordinate in ([5]) expressed with the help of x and y, i.e. formally we take f(x,y) — z(£,tp). 
According to our presumptions, this quantity should be negative for heptagonal defects and positive for pentagonal 
defects. 



A. Heptagonal defects 



Comparison of ([5]) and (PXT]) implicates 

*^_fM? =li (43 ) 

so 



f(x,y) = -\Jx 2 + y 2 - a? = -\Jr 2 - a 2 . (44) 
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After making up derivations d xx j ', d yy f and d xy f we get 



,,2 



K (7-2(1+7?) -fl 2 ) 2 ' (45) 

where we use the initial definition of 7?, given in (|51l) . As we clearly see, this expression is negative for arbitrary 
values of r. 



B. Pentagonal defects 



For the parametrization (1361) . we get 

f{x,y) 2 x 2 +y 2 



= 1, (46) 



so 

— \J x 2 + y 2 + a? = — 
a a 



f(x, y) = — a/x 2 + y 2 + a 2 = — \/r 2 + a 2 . (47) 

a 

Making up the corresponding derivations, we get using (|42|) 



,..2 



K (r 2 (l + ? /) + a 2 ) 2 ' (48) 



so the Gaussian curvature is all the time positive for the pentagonal defects. 



V. SOLUTION OF THE DIRAC EQUATION 



The solution of (|T9]) . (l20l) for heptagonal and pe ntag onal defects is introduced and the local density of states is 
calculated here. The linear elasticity theory [1J, [l6| is used. For the numerical calculations of LDoS, method 
described in [T3] is exploited. 



A. Heptagonal defects 



The form of (fT9 l) . (f20 l) will be 



— (j — $ cosh £) y tanh £ + r\u = E^/g^v, 



(49) 



- d^v - {j - $ cosh 2 £) J tanh 2 £ + r)V = E^/g^u, (50) 

where 

j = (j + 1/2- a v )/a, 5 = $/a, rj = c 2 /a 2 . (51) 

7] <C 1 is a dimensionless parameter which contains the elasticity properties of the initial graphene plane. Due to 
these properties, the defects can be interpreted as small perturbations in that graphene plane. In case of finite 
elasticity, we can do an approximation 77 ~ \fvt, where e < 0.1 [10| . By this way, the elasticity is described by a 
small parameter e. Its value is usually taken between 0.01 and 0.1. If we do some neccesary corrections of the gauge 
field lu^, for e — > we get the metric of the cone. 
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Suppose now E = 0. This energy corresponds to so-called zero-energy mode which is appropriate for the electron 
states at the Fermi level. Then, the solution of (l4"9l . (|5U|) is 



Uq 



Mf) - C '(A(0 + * ^0"-* ) ^ P (*^) , (.3, 



where 



k = t/T+v, A(0 = A/fc 2 cosh 2 £ - 1 (54) 



and C, C" are the normalization constants. 

For the nonzero values of E the solution can be written as in [llj 

5(0 = 5o(OW(0, «(0 = *to(0V(0, (55) 

where 

U(£) = W (0) (C) + eW (1) (0 + ' ' ' + e"W (n) (C) (56) 

and 

V(0=V( 0) (e) + £VW(0 + ---+e"VW(0, (57) 

e = pp 5 - and we take H = up = 1. Here n is an integer number and it is chosen according to our requirements on the 
precision. After substitution of this approximation into (|49|) and (|50p we get 

ddA = eA(£)V^, a e V = -eA(0W^. (58) 

u wo 

Putting W<°) = 1, V<°) = and i = 0, 1, n-1, in a numerical way, the form of the solution will be 

( 

U« +1 \0= /'v»(C)A(C)^dC, (59) 
o 

S 

V< i+1 >(0 = - / Z/W(C)A(C)^dC- (60) 
o 

The local density of states is for given £o defined as 

LDoS{E) = u 2 {E,^) + v 2 {E^ ). (61) 

For its evaluation we have to calculate the normalization constants C,C. They differ for different values of E. For 
unnormalized solutions 5'(£)>5'(£) of (|4"9")l . (f5C7|) and each value of energy, 

1/C 2 = 1/C' 2 = J (5'(0 2 +5'(0 2 )^. (62) 
o 

Since for ^ ma x — oo, the integral diverges, some finite value of ^ m ax in some interval, which is of particular interest, 
is needed. In this work, we take £ maa; = 2.5 and £ ma x = 2. It follows from the parametrization (1211) that for given £, 
the corresponding distance is r — acosh£, which means that for a — 1A we have r max = 6.13A, resp. r max = 3.76A. 
These values are of the same order as the size of the Brillouin zone which is formed by the single hexagons. Each 
atom in the hexagon lies at the distance 1.42 A from its nearest neighbours [TH, [l9j]. This is the main principle of the 
tight-binding approximation [20| in which we only admit the influence of the nearest neighbours. 
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FIG. 1: LDoS as a function of E € (—0.5,0.5) and £ £ (0,2.5) for 1-heptagon defects (left part) and 1-pentagon defects (right 
part) for B = 0; e = 0.01. 



B. Pentagonal defects 



The form of (fT9l) . (1201) is 



(i - $ sinri £) \f cotli £ + r}U = £\ 



(63) 



d^v - (j - $ sinh 2 £) v/ coth 2 £ 



77U 



In case of E — 0, the corresponding solution is 

u (f) = C(A(0 + fccoshO fc ^ + ^ 



sinh £ 



A(0 + cosh£ 



J / $A(0coshf 
ex P o 



(64) 



(65) 



where 



fc = + A (0 = V fc2 sinh2 f + L ( 67 ) 



For calculating the solution for nonzero values of E and the local density of states we use the same procedure as 
for the heptagonal defects. 



C. Local density of states 



In Figs. 1-3, LDoS as a function of energy and the parameter £ is presented for the surface with the defects formed 
by 1 polygon. In all these figures, we choose j = in (f5"Tj) and e = 0.01 in the expression for 77. We see the evidence 
that for growing B or £ m ax, the LDoS is decreasing and the decrease is faster for the pentagonal defects. If we took 
larger £ max , LDoS would go to zero with an exception of a small number of energies for which we would get plane 
waves. But larger values of £ are unphysical because of limited interval of validity of the tight-binding approximation. 



FIG. 2: LDoS as a function of E £ (—0.5,0.5) and £ £ (0,2.5) for 1-heptagon defects (left part) and 1-pentagon defects (right 
part) for B = 0.5; e = 0.01. 




FIG. 3: LDoS as a function of E £ (—0.5,0.5) and £ £ (0,2.5) for 1-heptagon defects (left part) and 1-pentagon defects (right 
part) for B = 1; e = 0.01. 



The values we chose enable us to compare LDoS for both kinds of defects for a small perturbation, where the 
difference between both approximations is not too large. Here we chose e = 0.01, but there aren't any significant 
changes for the LDoS if we let e to grow up to 0.1, as we can easily see from the graphs in Figs. 4 and 5, where we 
compare LDoS for different values of magnetic field at a particular distance for e = 0.01 and 0.1. 



VI. LANDAU STATES 



In H and [13, the Landau states for the researched defects are calculated for the conical and the hyperboloidal 
geometry for the positive curvature. In case of the conical geometry, where the corresponding parametrization is 

(r, tp) —¥ (r cos ip, r shop, cr) (68) 

(c is an arbitrary constant), we have, in agreement with the notation of |9j, 2 kinds of the Landau states. Let us first 
denote 

»=V±w i-444- < 69) 
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-0,5 0,0 0,5 -0,5 0,0 0,5 

E E 

FIG. 4: LDoS as a function of E £ (—0.5,0.5) for 1-heptagon defects (left part) and 1-pentagon defects (right part); various 
values of B are used, £ = 1.5, ^ ma x — 2; e = 0.01. 



TABLE I: Landau states for pentagonal defects for v < 0, B — 0.5 and j = — | 



n 







± 0.89, ± 1.18 


1 


± 1.34, ± 1.55 


2 


± 1.67, ± 1.84 


3 


± 1.95, ± 2.10 



Then the first kind coincides with the Landau states of the planar graphene [2l| . It corresponds to v > and the 
energy levels are 

E n = ±VWn, n = 0,l,2,..., (70) 

where B means the magnetic field. 

The second kind coincides with v < 0, the energy levels are 



E n = ±V2Bjn-v + -, n = 0,1,2,... (71) 

and the list of these Landau states for different values of the magnetic field is presented in Tables 1 and 2. Because 
j = — \ in (l69l) corresponds @ to j = in (fl"9|) , (|20|) . we are looking for the corresponding energy levels. 

Let us compare the Landau states for the hyperboloidal geometry with the values calculated for cone and graphene. 
For this purpose, we do an extension of the interval of energies for which we calculate the LDoS. The result we see 
in Fig. 6. We see that for pentagonal defects, the values ±1.67 and ±2.19 from Tables 1 and 2 can be found by this 
method. The reason for the presence or absence of other peaks could be the incompleteness of the list of Landau 
states for hyperboloidal geometry and a low magnitude of some peaks. It is also possible that some of Landau states 
characterizing the conical and planar geometry don't appear in case of the hyperboloidal geometry. For the negative 
curvature, the appropriate energy levels have similar position but are shifted to the left. 
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FIG. 5: LDoS as a function of E £ (—0.5,0.5) for 1-heptagon defects (left part) and 1-pentagon defects (right part); various 
values of B are used, £ = 1.5, (, m ax — 2; e = 0.1. 



TABLE II: Landau states for pentagonal defects for v < 0, B = 1 and j — — | 



n E. 






± 1.25, ± 1.66 


1 


± 1.89, ± 2.19 


2 


± 2.35, ± 2.59 


3 


± 2.75, ± 2.96 



Comparison of our results with the Landau states as expected for multilayer graphene [22j is also interesting. These 
results differ in the way that for zero-energy states (n = 0, see [22]), there are nonzero Landau states with low 
magnitude. For higher n, the Landau states are calculated with the help of the approximation formula (using the 
notation of (22J) 



( m ) ( mir 2 o mir \ 1 — (—1) N „ . mir n(n+l)Ti 



E K Z> , = - cos + cos — r 72 ± , / i — '—^ cos 4 — + — ^ — ' - , (72) 

"- M 1 A^ + l N + l 2{N+\)) 12 y N+l 12 2{N+1) 4 7l 2 sin 2 



where = \J^Bjo, N is the number of the graphene layers and 70 = 3, 71 = 0.4, 72 = —0.02. As we can see in 
Tables 3-5, the results acquired using this formula are higher in order then the Landau states calculated in this 
paper for the hyperboloid geometry. So, the Landau states for multilayer graphene and for the hyperboloid geometry 
are completely different. The reason may consist in considering the influence of atoms in neighbor and next- neighbor 
layers in the multilayer graphene which don't exist in case of the one-layer graphene. 



TABLE III: Landau states for 3-layer planar graphene, m = 2 



n 


B = 0.5 


E n 


B = 1 





± 16.87 




± 33.75 


1 


± 29.23 




± 58.46 


2 


± 41.33 




± 82.67 
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TABLE IV: Landau states for 4-layer planar graphene, m = 3 



n 




E n 






B = 0.5 




B = 1 





± 14.75 




± 29.50 


1 


± 25.55 




± 51.10 


2 


± 36.13 




± 72.26 


TABLE V: Landau states for 5-layer planar graphene, m = 4 


n 




E n 






B = 0.5 




B = 1 





± 13.78 




± 27.56 


1 


± 23.87 




± 47.73 


2 


± 33.75 




± 67.50 



VII. CONCLUSION 



Wc have studied the electronic structure of disclinated graphene in the vicinity of heptagonal and pentagonal defects 
depending on the kind of the curvature (negative or positive). Hyperboloidal parametrizations (|21)) .(|36 |) were assumed 
after rejection of the conical metric. The continuum field-theory gauge model was used in which the disclinations 
are incorporated using the vortex- like potential (p0|) . (|3Tj) for calculation the components of the metric. The arising 
fictious flux was compensated by the gauge flux of spin connection field (|15I) . ([T5| . The potential ([50)1 . (T5TT) also 
results in the dependence of the corresponding Dirac equation on the Frank index a which includes the number of 
defects. The defects are involved in (|5Tj) with help of the parameter e which comes from the elasticity properties of 
the graphene. 

Next we incorporated an uniform magnetic field (|34 |) . (|4T[) which can significantly influence LDoS. It is calculated 
from the solution of the Dirac equation which we get numerically with help of the extension of an analytical solution 
for zero-energy modes ([52")) . ([55)1 . (pjj) . . 

In all presented figures, the behaviour of the LDoS is compared for both kinds of defects. For very small values 
of e, this behaviour is for both kinds of defects similar, but it is more spread for heptagonal defects. As we found, 
the varying value of e does not affect LDoS significantly. It shows that with a little change of elasticity the change 
of LDoS is also very small. It could be done by big rigidity of graphene structure as well. We also compared the 
resulting Landau states with the theoretical prediction coming from the corresponding values for graphene plane and 
conical metric. We see that the different geometrical structure influences the position of the Landau states. 



'.0 




FIG. 6: Landau states for heptagonal (left part) and pentagonal defects (right part) for various values of B; (, ma x = 2, £ = 2, 
e = 0.01. 
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To conclude, the presented results could have a large potential use for next calculations of metallic properties of 
carbon nanostructures which serve as a wide spectra of electronic devices. Let us mention the significance of the 
zero-energy modes. Generally, they always appear as a solution for disclinated graphene in presence of the magnetic 
field [23| and they play a key role in explanation of anomalies, paramagnetism, high-temperature superconductance 
etc. 

We have to stress that we assumed defects in which appeared only 1 heptagon or 1 pentagon. For the higher number 
of polygons in defects, the calculation is more complicated especially for heptagons, because in contrast to pentagonal 
defects, problems with the geometrical interpretation occur. It is useful to do next calculations for more complicated 
forms of defects. 
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